Employing trapped cold ions to verify the quantum Jarzynski equality 



00 

o 
o 

<N 

< 

S 1 
Qh: 
+-> 
a 

3 . 
a*. 



> 

m 
m 
q 

od 
o 
oo 
o 



Gerhard Huber and Ferdinand Schmidt-Kaler 

Institut fur Quant en- Informationsverarbeitung, Universitdt Ulm, Albert- Einstein- Allee 11, D-89069 Ulm, Germany 

Sebastian Deffner and Eric Lutz 
Department of Physics, University of Augsburg, D-86135 Augsburg, Germany 

(Dated: August 4, 2008) 

We propose a scheme to investigate the nonequilibrium work distribution of a quantum particle 
under well controlled transformations of the external potential, exploiting the versatility of a single 
ion in a segmented linear Paul trap. We describe in detail how the motional quantum state of 
a single ion can be prepared, manipulated and finally read out to fully determine the free energy 
difference in both harmonic and anharmonic potentials. Uniquely to our system, we show how an 
ion may be immersed in an engineered laser-field reservoir. Trapped ions therefore represent an 
ideal tool for investigating the Jarzynski equality in open and closed quantum systems. 
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Nonequilibrium phenomena at the nanoscale are dom- 
inated by fluctuations and by quantum effects. The in- 
terplay of nonequilibrium thermodynamics and quantum 
theory is hence of fundamental interest. Only a decade 
ago Jarzynski published a major discovery in classical 
nonequilibrium physics [1], relating the free energy differ- 
ence AF after a given transformation to the probability 
distribution of the total work W done on the system: 



AF = —kT ln(e 



-W/kT\ 



(1) 



where {e~ w / kT ) = J dWe~ w/kT P(W) is the average ex- 
ponentiated work and k denotes the Boltzmann constant. 
This remarkable equality highlights the role of work fluc- 
tuations and provides a generic way of computing the free 
energy difference for any transformations, quasistatic or 
not, once the work distribution P(W) is known. Most im- 
portantly, the Jarzynski relation allows to determine AF 
even in the case of arbitrarily fast transformations, when 
irreversible thermodynamics is not applicable. Prior to 
the discovery of Eq. (pp), the determination of the free 
energy difference in such far from equilibrium conditions 
was not believed to be possible [2]. Recently, the classi- 
cal Jarzynski equality and its generalization by Crooks 
have been successfully tested in single-molecule stretch- 
ing experiments [4|, [5[. Later, the work distribution was 
recovered from repeated measurements of the mechanical 
work done on a colloidal particle in liquid solution [6] . 

The situation is much different at the quantum level. 
So far only studied theoretically, the Jarzynski equation 
holds in its classical form for closed quantum systems 
Q , while quantum corrections appear in the case of open 
systems due to the coupling to an external reservoir [8[. 
Further difficulties arise when considering the quantum- 
mechanical nature of work [9] and the question of how to 
measure it experimentally. It is evident that the classi- 
cal definition of work as force times displacement cannot 
simply be taken over unmodified. It has recently been es- 
tablished that work is actually not an observable in the 
usual sense, as it is not given by an expectation value 
of some Hermitian work operator, but rather by a time- 
ordered correlation function [10[. On the other hand, the 



problem of how to determine quantum work still remains 
unsolved, explaining the absence of an experimental veri- 
fication of the Jarzynski equality in the quantum domain. 

We show in this Letter how to experimentally measure 
nonequilibrium work using a single ion in a segmented 
linear Paul trap. A unique property of ion traps is the 
possibility to study the quantum Jarzynski equality, as 
well as the quantum Crooks relation for systems 
that are either isolated or coupled to tailored quantum 
environments using reservoir engineering [IH, Il3j ]. Single 
ions trapped in radiofrequency Paul traps are quantum 
nanosystems with remarkable properties. They can be 
laser cooled to very low temperatures, reaching to the 
motional ground state in the potential. Arbitrary quan- 
tum states can be prepar ed, manipulated and measured 
with high fidelity [14], LL5[ . Using the so-called electron 
shelving method, the quantum state is revealed with close 
to unity detection efficiency. The use of a segmented 
trap further allows for engineering a vast variety of time- 
dependent trapping potentials and hence the implemen- 
tation of different model Hamiltonians. In the follow- 
ing, we generalize the detection methods for the motional 
state pil 03 m order to realize an efficient filter for vi- 
brational number states. We show that trapped ions are 
not only good candidates for quantum computing, but 
may also allow to experimentally approach the emergent 
field of quantum thermodynamics. 

Quantum Jarzynski equality. We begin by consider- 
ing an isolated quantum system whose time-dependent 
Hamiltonian is varied from an initial value H(0) to a fi- 
nal value H(t). We denote by <\> l n and E l n the respective 
eigenfunctions and eigenvalues of the Hamiltonian H(t) 
at any given time t. We further assume that the system 
is initially thermalized at temperature T. The free en- 
ergy difference AF between final and initial state is then 
given by the Jarzynski equality, Eq. (pQ). The probability 
distribution of the random work W is given by [lOj], 



P(W) 



S[W-(E^-E° n )]P^ n P^ 



(2) 
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FIG. 1: Transition probabilities Pm, n f° r quantum numbers 
n = 2, m = to 8 as a function of the inverse switching 
time t -1 for a linear transformation of oo 2 (t) with oo(0) = 
0.1 MHz and uj(t) = 3.0 MHz. The inset shows the transition 
probabilities Pm,2 f° r a transformation time of r —1 jjls. 



where P° = (1/Z ) exp(-£°//cT) is the initial (ther- 
mal) occupation probability and P^ n are the transition 
probabilities between initial and final states n and ra, 



pr 



dx / d^^ n *(x)U(x,xo',T)^n( x o) 



(3) 



Here U(x,xq]t) is the propagator of the quantum sys- 
tem. The physical meaning of Eq. ([2]) is clear: the total 
work done during a given transformation of the Hamil- 
tonian is obtained from the energy difference between 
final and initial eigenstates, — E®, averaged over all 
possible initial and final states. Equation ([2|) shows in 
addition that the randomness of the work stems from 
the initial thermal distribution P®, and from the quan- 
tum nature of transitions between states, as described by 
Pm n . The origin of work fluctuations is therefore of both 
thermodynamical and quantum-mechanical nature. The 
free energy difference AF can be evaluated for an arbi- 
trary transformation of the Hamiltonian with the help of 
the Jarzynski relation, once the work probability density 
P(W) has been determined. We next describe a method 
to realize a quantum nonequilibrium situation for a sin- 
gle trapped ion in a linear Paul trap and how to measure 
its corresponding work distribution. 

Harmonic ion trap. Linear Paul traps are character- 
ized by a strong dynamical confinement in the radial di- 
rection (yz plane) and electrostatically bound in the ax- 
ial direction (x axis). With a radial confinement much 
stronger than the axial, we will restrict ourselves to the 
axial external degree of freedom. Near the center of the 
axial potential, the confinement is harmonic and the axial 
frequency uj can be varied in time by changing the con- 
trol voltages [18]. The quantum state of motion along 
the axial direction can be described by the Hamiltonian 




H(t) = 



P 

2M 



M 



-L0 2 (t)x 2 



(4) 



Work / tia^ 

FIG. 2: Work probability distribution, Eq. fl2}, for a linear 
change of the harmonic potential uj 2 {t) with u;(0) = 1 MHz 
and uj(t) = 3 MHz, for n = 1. Shown are the adiabatic case 
r — > oo (□) and two fast transformations with r = 0.1 {is 
(■) and r = 0.05 /is (•). For the latter, deviations from 
the adiabatic case are clearly visible. Even negative work 
processes arise: the probability contribution at W = —Hujo 
mainly stems from the transition n — 2 — > m — 0, which 
occurs with probability P^o = 10% (for r — 0.05 fis). This 
contribution can be readily tested by the number state filter. 



where M is the mass of the ion. For this simple potential, 
the nonequilibrium work distribution ([2j) can be studied 
analytically [19|. Besides the external, motional degree 
of freedom, the ion provides an internal, electronic level 
scheme. In our case, we consider a A-system comprising 
the ground state S1/2 and two excited states P1/2 and 
D 5 / 2 . The Pi/2 state rapidly decays into the S1/2, thus 
providing a high spontaneous photon scatter rate used 
for fluorescence detection. Laser induced transitions from 
the ground to the metastable D 5 / 2 state are induced on 
the narrow quadrupole transition (linewidth <C w{t)), 
if the spectral bandwidth of the Si/ 2 -D 5 / 2 exciting light 
field is small compared to the sideband structure. Co- 
herent laser pulses on this narrow-band optical transition 
allow to exploit and to store the motional quantum state 
information in the internal quantum states. 

The experimental measurement protocol of the work 
distribution consists of four consecutive steps: 

(I) The ion is first prepared in a thermal state with 
mean phonon number, n = (exp(fkJo / kT) — in the 
electronic ground state S-level by laser cooling and opti- 
cal pumping. We prepare this state deterministically by 
resolved-sideband laser cooling [2(j into the vibrational 
ground state |n = 0) and subsequently allowing the ion 
to heat up for a certain time without laser cooling. As 
the heating rate of the ion within the trap can be pre- 
cisely measured, this procedure is favorable for very low 
values of n. An alternative method, suited for higher 
values of n, is Doppler cooling on the S1/2 to P±/2 transi- 
tion. Varying the detuning of the cooling laser from the 
atomic resonance results in different thermal states with 
mean phonon numbers down to the Doppler limit. 

(II) In the second step, we measure the initial phonon 
number n using the filtering scheme described in detail 




FIG. 3: Levels and transitions involved in the filtering scheme. 
Shown are transitions driven by the initial 7r-pulse (i) and the 
successive 27r-red sideband (ii) and 27r-blue sideband pulses 
(iii), respectively. For |ratest), these pulses induce perfect 
7r and 27r-transitions between the metastable D 5 / 2 and the 
ground state level S1/2 (left side). No fluorescence is ob- 
served when the ion is exposed to resonant light on the S 
to P transition. For any other |n), the transfer will be im- 
perfect (dashed, right side) and there will be population in 
the S. Thus, the excitation to the P level is successful and we 
observe the emission of fluorescence photons. 



below. In this way, we determine the initial energy eigen- 
state E® (from spectroscopy measurements). 

(III) In the third step, we transform the trap poten- 
tial from an initial value uj(0) to a final value 00 (r) This 
changing potential will in general modify the ion's mo- 
tional state into a nonequilibrium state, while its inter- 
nal, electronic state remains unaffected. For simplicity, 
we consider a linear variation of the confining axial po- 
tential uo 2 (t) from a; 2 (0) to oo 2 (r). Figures [T] and [2] show 
a numerical evaluation, based on the results of Ref. [191 ] . 
of the transition probabilities (|3]) and the work distribu- 
tion ([2|) for realistic experimental parameters and differ- 
ent transformation times r. 

(IV) In the last step, we measure the new phonon num- 
ber m using the filtering scheme and determine the final 
energy eigenstate E^. The distribution of the nonequi- 
librium work, W = E^ — Eq. ([2]), is then recon- 
structed by repeating the measurement sequence. By 
evaluating Eq. (pQ) for adiabatic and nonadiabatic pro- 
cesses, we can verify the Jarzynski equality. 

Filtering scheme. A sequence of laser pulses on the 
narrow S to D transition is applied to the ion, coherently 
processing its internal and external degrees of freedom 
[2l| . We tailor this pulse sequence such that the ion will 
end in the metastable D 5 / 2 state with certainty if the vi- 
brational quantum state was |m te st)- Subsequently, the 
ion is illuminated with light resonant to the S1/2 to P1/2 
transition. If we observe no fluorescence, the ion is mea- 
sured in the D state. However, for vibrational states 
different from | mtest)? the laser pulse sequence leads to 
a superposition state, a\S) + /3\D), such that there re- 
mains a non- vanishing probability \/3\ 2 of projecting the 
superposition into \D), and thus observing no fluores- 
cence. We therefore repeat the procedure a few times 
such that a high quality of the filtering procedure is en- 
sured. Considering the evolution of the quantum state 
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Phonon number 

FIG. 4: Probability for detection of zero-fluorescence, called 
filter transmission, after N = 1, 2, ... 10 pulse/detection trial 
cycles. With increasing number of trials (from back to front) 
the Fock state to be tested for, here mtest = 3, always remains 
dark under excitation with resonant light on the S1/2 to Pi/ 2 
transition. All other states show an exponential decrease of 
zero-fluorescence detection probability. 



itself, the influence of the above sequence reminds of the 
operating principle of a filter: its projective 'transmis- 
sion' is unity for a certain input state | mtest) an d zero 
otherwise. We adapt the laser pulse sequence timing to 
reach all relevant eigenstates \n) and |m) with the filter. 

The crucial requirement for a well-suited filter proce- 
dure is to ensure the non-zero fluorescence outcome for 
all states but |ra te st)- It is sufficient to design the number 
state filter to have high suppression factors in a vicinity 
of I m t est) 7 since n is rapidly vanishing for high values 
of |m — n\. Varying the duration of the pulses, we use 
the following procedure, see Fig. [3j Starting from state 
\S, mtest), we apply a 7r-pulse on the first red sideband 
leaving the ion in \D, mtest — 1)- As the Rabi frequency 
^n,n ; between vibrational states n and n r depends on 
both initial n and final n' , the laser pulse does not com- 
pletely transfer other vibrational states to the D state. 
When we expose the ion to resonant light on the S to P 
transition, zero-fluorescence is observed if the ion was in 
latest), but for other vibrational states there is a certain 
probability to observe fluorescence. To sharpen the dis- 
crimination, we apply 27r-pulses on the red and blue side- 
band, respectively, interleaved by a fluorescence detection 
trial after each pulse. Again, the 27r-pulses and detec- 
tions leave the ion in the dark state rarest — 1)? but 
yield a non-zero fluorescence signal for all other states. 
This resembles so-called trapping states which have been 
investigated in cavity QED experiments [22[. This prob- 
ability for zero-fluorescence detection decreases exponen- 
tially with the number TV of pulse/detection runs, as 
shown in Fig. [H for a wide range of states mtest > 3. 
After N = 10 runs, all probabilities but for |m te st) drop 
below 5 %. As no coherence is remaining after each de- 
tection interval, the scheme has modest requirements on 
the phase stability. The driving laser needs preserving 
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phase only for one single 27r-pulse, but not during the 
entire filter sequence. The two lowest vibrational states 
are treated even simpler: For m te st=0, a 7r-pulse on the 
carrier transition brings the ion into \D,0). Successive 
red sideband 7r-pulses do not affect this state, but fluores- 
cence is observed with non-zero probability for all other 
states. This scheme has been proposed for stochastic 
cooling pilGj]]- For mt es t=l, the carrier pulse is simply 
replaced by a red sideband 7r-pulse |S, 1) — > |D,0). Then 
the procedure continues as for m t est=0. The length of 
the pulses is specific for each choice of m te st 7 varying the 
pulse allows to access measurements over a wide range of 
vibrational levels. 

To estimate the time taken by one experimental cycle 
from preparation to detection, we assume a few 10 [is 
for sideband pulses and a few hundred /is fluorescence 
detection time; one cycle with multiple filtering iterations 
will then take less than 10 ms, short compared to the 
lifetime of the D-state (1.2 s for 40 Ca + ). To also assure 
that unwanted dissipative effects do not introduce errors, 
the trap's heating rate needs to be smaller than 1 phonon 
within the cycle time. Traps with much lower rates have 
been reported. The statistical error of the values P^ n n is 
further reduced by repeating the measurement sequence. 

Designing the bath properties and the potential shape. 
As discussed before, single trapped ions are highly iso- 
lated from external reservoirs. However, it has been 
shown theoretically [12[, and also in first experiments 
[l3^ that it is possible to introduce a coupling between 
the ion and an artificially laser-induced bath. The vari- 
ation of laser frequencies and intensities allows one to 
engineer the coupling and select the master equation de- 
scribing the motion of the ion. Here, the interaction is 
mediated by sideband transitions between the Si/ 2 and 
D 5 / 2 level, see Fig. 03 For example, a zero-temperature 
reservoir can be implemented by a light field tuned to the 



(cooling) red sideband transition. A large variety of other 
tailored reservoirs, such as squeezed baths, can be gener- 
ated as well. Within the framework of this proposal, it is 
therefore possible to investigate nonequilibrium transfor- 
mations of open systems with tailored couplings between 
system and reservoir. In particular, the distribution of 
the heat exchanged with a reservoir can be determined 
using the same measurement protocol by keeping the trap 
frequency constant, that is, performing no work. 

Exploiting the flexibility provided by a segmented trap 
design, it is moreover feasible to investigate anharmonic 
trapping potentials. Especially in situation of a nonadi- 
abatic transport along the segments of the trap [18| , the 
ion is shifted out of the harmonic center of the electric 
potential and experiences nonharmonic potential contri- 
butions [2O]. For future work, one might include forces 
by laser light on the ion, which depend on its internal 
electronic state, investigating the influence of quantum 
thermodynamics on qubit gate operations [231 ]. 

In conclusion, we have shown how the quantum 
Jarzynski equality can be experimentally investigated us- 
ing a single ion in the time-varying electrical potential 
of a Paul trap, for both open and closed quantum sys- 
tems. Our proposal is based on the state-of-the-art in 
many laboratories working with single trapped ions and 
uses realistic parameters. Experiments with such a de- 
vice would allow to shine more light on the amazing in- 
terplay of quantum mechanics and thermodynamics. 
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